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We examine spectral properties of doped holes dressed with surrounding spin cloud in the t- 
J model. These composite-hole excitations well characterize prominent band structures in the 
angle-resolved photoemission spectrum. In one-dimensional (ID) case at half-filling, we identify the 
composite operators that separately pick up the spinon and holon branches, respectively. After hole 
doping, we find that the composite hole excitations with string-like spins tend to be localized at 
k = 7t/2 in the momentum space. This means that such composite excitations should be actual 
electronic excitations, since the spinon and holon branches merge together at this momentum. In 
2D case, we find that the composite excitations with more non-local spin fluctuation have stronger 
intensity near the Fermi level. The composite band structure along diagonal (0, 0)-(7r,7r) direction 
in 2D has some similarity to that in ID, and such non-local spin fluctuation plays an important role 
on the formation of the pseudogap in high-T c cuprates. 

PACS numbers: 74.25.Jb, 71.10.Fd, 74.72.Kf, 74.25.Ha 


I. INTRODUCTION 

The t-J model is one of prototypical examples for un¬ 
derstanding strongly correlated electron systems. The 
interplay between spin and charge degrees of freedom is a 
key ingredient of this model. In spatially one-dimensional 
(ID) case, this model (or equivalently the large-U Hub¬ 
bard model) has been fully understood by the exact and 
numerical approaches [l,|2]. Therein, the electronic band 
dispersion splits into the spinon and holon branches, and 
this is called the spin-charge separation. They are col¬ 
lective excitations, not the quasiparticles. In 2D case, it 
is believed that this model captures anomalous features 
of hole-doped high-T c superconductivity in cuprates Q. 
Although it is difficult to solve it exactly, various types 
of numerical studies are proceeding. The value of this 
model as the prototypical and minimal correlated sys¬ 
tem is still not lost, even though roughly 30 years have 
passed already after the discovery of high-T c cuprates. 

Generally speaking, the spin and charge degrees of free¬ 
dom are separated in ID, while they are interacting with 
each other in 2D. However, this statement is sometimes 
misleading, and for example in the context of the slave 
boson theory the phase string connects the spinon with 
the holon excitation Q. The notion of the spin-charge 
separation and coupling may depend on what we are go¬ 
ing to observe and on how to treat the spin and charge 
degrees of freedom. In this paper, we would like to revisit 
the coupling of spin and charge degrees of freedom in ID. 

Mathematically, the fundamental operator of the t-J 
model corresponds to a fermionic particle with no double 
occupation at each site. This operator does not obey the 
exact fermionic anti-commutation relation. Owing to this 
constraint, the equation of motion for this excitation nec¬ 
essarily induces more non-local spin excitations [Ej], and 
the equations of motion do not close. A doped hole into 
the Mott insulator becomes a composite particle dressed 
with the spin excitation. This mathematical structure is 
common to both of ID and 2D cases, and thus it is nat¬ 


ural to consider that this non-locality is related to the 
phase string in ID and the spin cloud in 2D. The pur¬ 
pose of this paper is to examine the spectral properties 
of such composite excitation modes. Based on the com¬ 
parison between ID and 2D cases, we would like to also 
discuss about the anomalous band structure of high-T c 
cuprates. 

By this approach, it is possible to do more advanced 
spectroscopic calculation rather than the standard angle- 
resolved photoemission spectroscopy (ARPES). In terms 
of composite operators, the complicated band structures 
of correlated systems are composed of various types of 
composite particles. Different branches are represented 
by different composite operators, respectively. Even 
if the single-particle spectrum is heavily damped by 
the strong self-energy correction at a particular energy- 
momentum region, we may find the composite excitation 
that is well defined and has a long life time at the re¬ 
gion. In other words, we would like to find the best rep¬ 
resentation of the collective mode at given energy and 
momentum. 

The organization of this paper is as follows. In the next 
section, we introduce the technical details. In particular, 
we explain the physical meaning of composite excitations. 
In Secs. Ill and IV, we present numerical results for ID 
and 2D cases. Based on these results, we compare ID 
with 2D results, and finally we summarize our paper. 


II. METHOD 

A. Model and Fundamental Algebra 

We start with the t-J model. The Hamiltonian is de¬ 
fined by 

H = — tij Hi • Jij , 

i,j <ij> 


( 1 ) 
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where = (£j^,£j^) is the spinor representation of 
singly-occupied particle at site i, ^ ?cr = Ci, a ( 1 — n^-o) 
with spin cr, = E<t 4 ,<t 6 >> n i = 4 6 = n *,0 = 

4 (T o 6 , and Hi = 4<?6 with the identity matrix <7o and 
the Pauli matrices 3 = (oi, cr 2 , 03 ). This model includes 
electron hopping, and the superexchange interaction 
J. Here, the hopping term between the nearest neighbor 
pair of sites i and j is represented as 

t{j — tOLij . ( 2 ) 
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( 10 ) 


The Fourier transform of otij is given by 

a(k) = 2 cos/c, (3) 

for the ID chain, and is also given by 

a(k) = 2 (cos k x + cos fc^), (4) 

for the 2D square lattice. 

It is noted that the excitation ^ is not a simple fermion, 
and satisfy the following algebrae 


{6,4} = 1 + 

( 5 ) 

= -rii + a -Hi, 

(6) 


( 7 ) 


(8) 


where 1 = gq. As we can easily imagine, the first equality 
necessarily produces spin and charge fluctuation. As we 
go to higher order hierarchy of the Heisenberg equation 
of motion for an excitation operator, the fundamental 
excitation is dressed with such quantum fluctuation. 


B. Hierarchy of Composite Excitations and Their 
Propagators 


These excitations have large spectral intensities that are 
defined as the diagonal parts of (0| {T, 4^} |0) @. Thus 
if we look at the spectra of these composite excitations, 
we may say the origin of each prominent band structure 
in the ARPES spectrum. We have confirmed that the 
excitation originating from the exchange term, 4>2 = H a • 
3%, show similar behavior to T 2 , and hereafter we do not 
consider it. 

In this paper, we calculate the diagonal part of the 
following propagator matrix 


G °^ = < 0| * ^ + a- g + A '°> 


(0| T 


1 


t_ 

k,(7 0 J E 0 + H + ie * k ’ a 


%J0) 


( 11 ) 


where |0) is the ground state, and is the Fourier 
component of The spectra for the composite par¬ 

ticles are given by the imaginary part of the retarded 
propagator matrix 

A(k,w) = — lim — ImG R (k, cj). ( 12 ) 

e —>-0 7 T 

For instance, the usual single-particle spectrum is given 
by Hereafter, we particulary focus on 

(k) lo') j and A^++^++t (A), cj). 


We consider the Heisenberg equation of motion for 
£i(t) = e~ lHt ^ie lHt (we take h = 1). This equation of 
motion does not close, and induces spatially more ex¬ 
tended composite excitations. The equation for ^(t) is 
given by 

* = ~t(l+ £f(t) 

+ (9) 

where we define = JA otijij- Then, we find the new 
excitations and in which a doped hole com¬ 

bines with the nearest neighbor spin fluctuation. When 
we next consider the equation of motion for these com¬ 
posite states, we find Hi • 3^f a and Hi • 3 (Hi • 3^) a . For 
instance, the important excitations up to L = 5 for 2D 


C. Numerical Conditions 

In ID (2D) case, we take 16 (4 x 4) lattice under the 
periodic boundary condition. In doped cases, we take 
2 holes (14 electrons), and thus the hole doping rate is 
S = 0.125. The Hamiltonian matrix is exactly diagonal¬ 
ized by the Lanczos algorithm for the calculation of the 
propagators for composite excitations. The propagators 
are also calculated by the kernel polynomial expansion, 
and we have confirmed the relevance between these meth¬ 
ods. The broadening factor is taken to be e = 0 . 2 1. Since 
the momentum resolution of each spatial direction is 7r/2 
owing to severe finite-size effect, we naturally interpolate 
the numerical data that the mesh size becomes 7 r/ 20 . We 
should note that this procedure does not produce new 
fine structures originating with lower energy physics. We 
just introduce this method so that our data are visible. 
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III. NUMERICAL RESULTS I: ID CASE 
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FIG. 1: ARPES spectrum (spectrum for the bare fermionic 
particle = £) 




FIG. 2: Angle-resolved spectrum for ^2 = n • 

Let us first examine the half-filled case. We present 
the angle-resolved spectrum for each component of 
in Figs. CDS The figure [ 1 ] for = £ corresponds to 
the standard ARPES result. Our numerical result of the 
ARPES spectrum is consistent with the past analytical 
and numerical results, and we clearly observe the spinon 
and holon branches with different energy scales J and 
£, respectively. Based on this consistency, let us par¬ 
ticularly focus on the spectra for ^ 2 , ^ 3 ", and in 

the following sentences. The doped holes in these modes 
strongly couple with various sizes of the spin strings. 
Thus, they may not be simple single-particle-type exci¬ 
tations, although they are still fermionic excitations. We 
expect that they would play crucial roles on the forma- 
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FIG. 5: Angle-resolved spectrum for 
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FIG. 6: Angle-resolved spectrum for T 4 + 



FIG. 7: Angle-resolved spectrum for 4^ 



FIG. 8: Angle-resolved spectrum for 4'J“ + 



® 

® 


® 

® 


<D 

<D 


® ®f® ® ® ® ® ® 

® ®(6® ®A® ® 


FIG. 9: Photoemission processes from AF-like background: 
(a) the electron removal process by 4^ and (b) the electron 
removal process by 4 r J + . 


tion of collective excitations such as spinon and holon 
excitations. 

At first, we examine the spectrum for 4^ = n • <r£ a 
in Fig. [2j We find that the main band exactly traces 
the spinon branch. This is quite unique, since such a 
property enables us to do more advanced spectroscopy 
in the sense that we can select one structure of compli¬ 
cated bands in correlated electron systems. The reason 
for such separable spectroscopy is quite simple. When we 
remove one electron from the AF insulating ground state, 
the domain boundary (spinon) of the AF state is created 
at around the position where the electron is removed. 
Since the electron removal by the operator T 2 counts the 
change in the near-neighbor spin structure as well as the 
created hole, the operator necessarily presents the spinon 
creation [see Fig. Ha)]. Therefore, we can regard T 2 as 
’the spinon operator\ Note that this is different from 
the spinon operator / in the context of the slave boson 
representation £ = ftf. 

On the other hand in Figs. 0] and [ 8 j the main bands 
for T 3 " and 4 / 4 trace the holon branch, and the spec¬ 
tral weight of the spinon branch is not strong. Thus, they 
are basically ’the holon operators\ It is noted that the 
covering of the holon branch by these types of operators 
becomes much better as the spin string becomes longer. 
The readers may wonder why they behave like the holon 
excitation, although they are dressed with many spins. 
The answer may originate in a fact that the spin string 
attached to ^3 and 4/^ + can regulate the spin struc¬ 
ture after the photoemission so that the resulting final 
state has just one hole. Let us consider the situation in 
which they are operated to the AF-like state, but the 
state contains one spinon bases owing to quantum fluc¬ 
tuation. Then, the operation of 4 / 3 " and 4f^ + to the state 
can remove the spinon as well as creation of the hole [see 

























Fig. lib)]. 

These results suggest that each branch in the ARPES 
data originates from a particular composite excitation. 
Interestingly, these spectral intensities are much stronger 
than that of the single-particle spectrum, and there is 
only weak continuum part. This weak continuum fea¬ 
ture really represents that these composite excitations 
are more important than the simple £ state. 

In comparison with T 2 , and the spectra 

for T 4 - , and as shown in Figs. [3j 

O EJ and 0 are more scattered, though a considerable 
amount of spectral intensity is still on the holon branch 
at 7 r /2 < k < 7 T. For these excitations, the multiple spin 
excitations are away from a created hole by the photoe¬ 
mission, thus they seem to be somehow different from 
T 2 , and 

B. Hole doping 
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FIG. 10: ARPES spectrum (spectrum for the bare fermionic 
particle Ti) 

Next we examine the doped case in Figs. [TOlfITl In 
this case also, the important composite excitations near 
the Fermi level are T 2 , Tg - , and We find that the 

spectral intensity concentrates on k = n/2 and uj/t = 0 
as the doped hole is dressed with longer spin string. At 
this momentum, the spinon and holon branches merge 
together. Thus, the angle-resolved spectra for T 3 " and 
have enough information for real electronic excita¬ 
tion as well as the phase string that connects spinon with 
holon. 

In the standard ARPES spectrum (Fig. fTO]) . we observe 
the metallic branch showing the Tomonaga-Luttinger liq¬ 
uid behavior. The spectral weight above the Fermi level 
originates in the holon branch at k > 7 r/ 2 . However, for 
the composite excitations T 2 , T 3 ", and T 4 + , the spectral 
intensity of this branch tends to be vanishing above the 
Fermi level. In particular, the spectral intensity for 
and is strong only at the Fermi level. On the other 
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FIG. 11: Angle-resolved spectrum for = n • 
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FIG. 12: Angle-resolved spectrum for T 3 = ft • cr£ ) cxa 



FIG. 13: Angle-resolved spectrum for = n • <7 (n • 
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FIG. 14: Angle-resolved spectrum for T 4 
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FIG. 15: Angle-resolved spectrum for T 4 + 
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FIG. 16: Angle-resolved spectrum for Tj 
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FIG. 17: Angle-resolved spectrum for 

hand, the spectrum for has considerable weight above 
the Fermi level. Therefore, the composite excitations be¬ 
have quite differently from the fundamental excitation 
Ti = £. Although we observe single metallic band in the 
ARPES spectrum, these results suggest that the charac¬ 
ter of the wave function changes at the Fermi level. This 
is because for instance T 2 corresponds to the self-energy 
correction to Ti in the standard field-theoretical termi¬ 
nology, and this means that the damping mechanism of 
the fundamental excitation strongly depends on the mo¬ 
mentum. We will later discuss about this feature again 
in comparison with the presence of the pseudogap in 2D 
cases and high-T c cuprates. 

The spectra for T 4 - , T 4 + , and seem to con¬ 
tain different information. Actually, we observe relatively 
strong spectral intensity at 0 < k < tt and u/t ~ 1. This 
shadow band seems to naturally connect to the holon 
branch at 7 r/ 2 <&< 7 r, and intersects with the metal¬ 
lic band at around k = 7 r/ 2 . Owing to this intersection, 
the band seems to split into two substructures. Later, 
we will also discuss more about this intersection in con¬ 
nection with the presence of the pseudogap in high-T c 
cuprates. 

IV. NUMERICAL RESULTS II: 2D CASE 

Let us move to discuss the band structure in 2 D case. 
We examine the hole doped case (6 = 0.125) shown in 
Figs. I22ll25l For comparison, we also show undoped case 
in Figs. UHHHJ We should be careful for the fact that the 
numerical data may involve strong finite-size effects. The 
excitations T 4 ^ wrap the overall lattice sites, and their 
motions are highly restricted. Thus, we only show the 
spectral data for Ti, T 2 , T 3 , and ^ 3 “. 

We examine the spectrum for Ti = £. The data are 
shown in Fig.[22j In the figure, the Fermi level is located 
at the origin of the energy cu = 0. Comparing this with 
Refs. [7, 8 ], we see that the result is reasonable. The non- 
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FIG. 18: ARPES spectrum (spectrum for 'hi = £) at half- FIG. 21: Spectrum for 'hj = ft • a (ft • <J^ Q! ) a at half-filling 
filling 
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FIG. 24: Spectrum for T 3 = n • for the two-hole state 
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FIG. 25: Spectrum for = ft • <7 (n • <r£ a ) a for the two-hole 
state 


interacting band is highly renormalized, and the band 
width becomes roughly half. Furthermore, when we look 
at the spectrum near the Fermi level, we find that the 
weight at k = (7t/2, 7t/2) and (7r, 0) is very strong. When 
we look at the branch from (0, 0) to ( 7 r, 7 r) through ( 7 r, 0), 
we find that the weight just above the Fermi level is once 
reduced at around ( 7 r, 0). This tendency seems to be con¬ 
sistent with the presence of the Fermi arc. We have also 
confirmed that the gap clearly opens at this momentum 
by introducing long-range hopping t' and t". Accord¬ 
ing to Ref. 0, the pseudogap should be also seen just 
above the Fermi level at around (7t/2,7t/2) (This gap is 
not seen experimentally, since the gap is located at the 
electron addition side). The branch from (0,0) to (77,7r) 
through (7t/2,7t/2) suffers from the presence of the gap 
owing to our finite lattice system, and thus unfortunately 
we cannot determine whether the gap exists. It is also 
noted that there is a considerable amount of broad spec¬ 
tral weight at around k = (77,7r) and —4 < u/t < —2. 


We also examine the spectra for T 2 , T 3 , and in 
Figs. 1231 (Ml and 1251 respectively. We find that T 2 and 
Tjj" have strong spectral intensity near the Fermi level, 
while 4/^ does not have the intensity near the Fermi level. 
Thus, the strong coupling between the doped hole and 
the spin cloud is crucial for the intensity even in 2D cases. 
Furthermore, the intensity near the Fermi level increases 
as the doped hole is dressed with larger and larger spin 
cloud. All of the excitations do not have the remarkable 
spectral intensity above the Fermi level. This feature is 
also roughly consistent with ID case, and strongly sug¬ 
gests that the character of the wave function changes at 
the Fermi level. 


V. COMPARISON BETWEEN ID AND 2D 
RESULTS 

It would be quite meaningful to compare the ID results 
with the 2D results along (0, 0)-(7r, 7t) direction through 
k = (7t/2,7t/2) 0. For instance, when we look at the 
spectra for 4 / 2 in ID and 2D cases, we see that they 
are quite similar with each other. This observation sug¬ 
gests that the large spectral intensity at k ~ ( 77 , n) and 
uj ~ —4 1 in 2D originates in the spinon-like excitation 
observed in ID. Originally, if we do not consider the 
spin fluctuation, there is no spectral weight at around 
k = (7r,7r) and uj/t ~ —4. Thus, this is quite reasonable. 
However, if we look at ^3 and the intensity at 

k ~ ( 77 , 77 ) and oo/t ~ — 4 in 2D is much enhanced than 
that in ID. This may be due to dimensionality depen¬ 
dence on the coupling strength between spin and charge 
degrees of freedom. Since the ID metallic band originates 
in the holon branch before doping, in 2D case also the 
weight at around —4 < u/t < — 2 may be pushed above 
the Fermi level by the hole doping. This scenario is ac¬ 
tually consistent with our conclusion that the character 
of the wave function changes at the Fermi level. 

If the dimensionality dependence is less pronounced 
on the spectra of composite particles, we may say that 
the presence of the pseudogap in the 2 D t-J model and 
high-T c cuprates can be understood by considering the 
vanishing spectral weights of the composite excitation 
above the Fermi level in ID cases. I believe that this 
observation is a strong indication that the pseudogap is 
different from the superconducting gap. However, this is 
still highly hypothesized conjecture, and thus we need to 
examine more about the vanishing spectral weight above 
the Fermi level in ID and 2D in more physical stand¬ 
points. 


VI. SUMMARY 

Summarizing, we have examined spectral properties 
of composite excitations in the ID and 2D t-J models. 
The most important goal of this numerical work is that 
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these composite hole excitations well characterize promi¬ 
nent band structures in the angle-resolved photoemission 
spectrum. Therefore, it is possible to do more advanced 
spectroscopy rather than the standard ARPES spectrum. 
In ID case, we have identified that the spinon and holon 
operators are represented as T 2 , Tj, and After 

doping to the Mott insulator, the spectral weight of 
and concentrates on the Fermi level, and they are 
physical excitations. We have also found that the ID 
composite data are quite similar to those in 2D except 
for relatively large spectral intensity of 2D at k ~ (n, tt) 
and l u/t ~ —4. In both of ID and 2D cases, we have 


found that the band dispersions above and below the 
Fermi level have different characters, respectively. We 
have proposed that this feature would provide us quite 
important information on the formation of the pseudogap 
in high-T c cuprates. 
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